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1 Introduction 

A major research focus in economics and econophysics is on the distribution of wealth in societies 
at different stages of development. Wealth includes money, material goods and assets of different 
kinds. Knowledge of the monetary equivalent of the latter two components is required in order 
to quantify wealth. A related and easier to measure distribution is that of income. The major 
motivation of theoretical models is to provide insight on the microscopic origins of income/ wealth 
distributions. Such distributions are expected to provide good fits to the empirical data. In the 
context of incomes, Champernowne 2 has commented "The forces determining the distribution of 
incomes in any community are so varied and complex and interact and fluctuate so continuously, 
that any theoretical model must either be unrealistically simplified or hopelessly complicated." 
The statement highlights the desirability of finding a middle ground between the unrealistically 
simple and the hopelessly complicated. 

A number of distribution functions has been proposed so far to describe income and wealth 
distributions. Theoretical models based on stochastic processes, have been formulated to explain 
the origins of some of the distributions^"^. One proposed distribution, mention of which is found 
in economic literature, is the beta distribution^"^ . In this paper, we describe a simple stochastic 
model of wealth distribution and show that the beta distribution is obtained in the non-equilibrium 
steady state. 

2 Stochastic model 

In the model, each economic agent (can be an individual, a family or a company) may be in two 
states: inactive (E) and active (E ). We determine the probability distribution of the wealth of 
an agent randomly selected from a population of agents. Let the agent possess wealth M at time 
t. Increase in the wealth of the agent can occur in two ways: at a steady rate and at random time 
intervals. In state E, the agent's wealth increases at rate b m and in state E , the rate is given 
by b rn + j m . In both E and E , the agent's wealth decreases at rate k m M. The decay rate is 
proportional to the current wealth with k m being the decay rate constant. Transitions between 
the states E and E occur at random time intervals. The rate of change of wealth is governed by 
the equation 

dM . 

—j- = JmZ + b rn - K m M = f(M, Z) (1) 

where z — 1(0) when the agent is in the state E (E). Let pj(M,t) (j = 0, 1) be the probability 
density function for wealth distribution when z — j. The rate of change of the probability density 
is given by 
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dpj{ ™' t] = -^[f(M,z) Pj (M,t)]+Y^[W kj p k (M,t) - W jkPj (M,t)\ (2) 

where W k j is the transition rate from state k to state j. The first term in Eq.(2) is the "transport" 
term representing the net flow of probability density and the second term represents the gain/loss 
in the probability density due to random transitions between the state j and the other accessible 
state. One can define the activation and deactivation rates, k a and kd respectively, to be k a = W n i, 
and k d = W w . From Eq. (2), 

= -gj^{(b m -k m M)p } + k d p 1 -k a p (3) 
dp d 

~~dt = ~dM^ m + bm ~ k mM) Pl } + k a p - k dP i (4) 

with p = po + pi . The minimum and the maximum amounts of wealth possessed by the agent are 
given by M min = b m /k m and M max = (b m +j m )/k m . Define to = M/M max , m mm = M mm /M max , 
i~\ = k a /k m and r 2 = kd/k m . In the steady state, dpa/dt = and dpi/dt = 0. The steady state 
solution turns out to be the beta distribution 



(rn - m m i n ) ri 1 (1 - m) r2 1 
B(n,r 2 ){l - m mm ) r 



= 571 — ~ w+r2 -i ( 5 ) 



The normalization constant i3(ri,r 2 ) is 

B (^^)= (1-m ■ )r 1+ r a -l dm ( 6 ) 

Jm m i n V l,tj min) 

In Eqs. (5) and (6), r± > 0, r 2 > and m m i n < to < 1. Let m m in = and n > 1 and r 2 > 1. In 
this case, 

r(ri)r(ra) 

B(ri ' r2)= i>7T^) (7) 

is the well-known beta function. The mean wealth m av and the variance m var are given by 

n rir 2 1 



j m var , , 

ri + r 2 (ri+r2) ri + r 2 + 1 



(8) 



The quantities depend on the ratios n and r 2 rather than on the individual values of fc 0j fcrf and 



3 Results and discussion 

Societies are traditionally divided into three classes: poor, middle and rich. Figs. 1(a), (b), and 
(c) show the p(m) versus m distributions in the three cases. One can obtain similar curves when 
mmin 7^ 0. The Gini coefficient G, a measure of wealth inequality, is expected to be small for each 
separate class. For example, G = 0.2 in the case of Fig. 1(a) describing wealth distribution for 
the poor class. The two-parameter beta distribution is flexible and can take a variety of shapes. 
The precision in fitting data is, however, limited in this case. 



McDonald and Xu 10 have proposed a five-parameter generalised beta distribution 

a ' 6 ' C ' P > q) = b^B(p, q )(l + c( m ^ (9) 
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(a) 



(b) 



(c) 



Figure 1: Probability density function p(m) as a function of m for (a) n = 2, r 2 = 8, (b) 

n = 10, r 2 = 10, (c) n = 10, r 2 = 2 

where < y a < 6 a and is zero otherwise. Also, < c < 1 and b, p, q > 0. B(p, g)represents the 
normalisation constant. The beta distribution (Eq. (5)) is a special case of GB(y; a,b,c,p,q) with 
fn-min = 0, c = 0, a = 1, and 6=1. 



Table 1. Generalised beta distribution (GB) and its special cases 



5 parameter 



GB 



c = 



c = 1 



4 parameter 



3 parameter 



q = l 
a = — 1 



2 parameter 



Pareto 



GB1 



a = 1 



GB2 




a = 1 



betal 



GG 




a = 1 



lognormal 




beta2 



Singh- 
Maddala 




Daeum 



gamma 



Weibull 



Many well-known distribution functions are limiting cases of the generalised beta distribution 
GB. Some examples are shown in Table 1. The betal distribution reduces to the beta 
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distribution (m min = 0) with 6=1. GG refers to the generalised gamma distribution. The 
special cases of GB(y), GB1 and GB2 have been shown to outperform other distributions in 
providing good quantitative fits to the income data from various countries and segments of 
society. The beta distribution, considered in the paper, is a special case of GBl. 

In this paper, we have provided a stochastic model of wealth distribution leading to the beta 
distribution in the non-equilibrium steady state. It will be of interest to formulate stochastic 
models of generalised beta distributions GB, GBl and GB2. An understanding of the microscopic 
origins of income/ wealth distributions may provide insight on the policies required to ensure that 
the benefits of economic growth reach all sections of society. 
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